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PnThe iron-based LaFeAsOi-xF^c recently discovered by Hosono's group is a fresh theoretical challenge as a new class of high- 
itemperature superconductors. Here we describe the electronic structure of the material and the mechanism of superconductivity. 
^ We start with constructing a tight-binding model in terms of the maximally localized Wannier orbitals from a first-principles 
^ electronic structure calculation, which has turned out to involve all the five Fe 3d bands. This is used to calculate the spin and 
d tharge susceptibilities with the random phase approximation. The spin susceptibility has peaks around k = (7r,0), (0, tt) arising 
I from a nesting across disconnected Fermi surface pockets. We have then plugged the susceptibilities into the linearised Eliashberg 
equation. For the doping concentration x = 0.1, we obtain an unconventional s-wave pairing, which is roughly an extended s in 
that the gap changes sign between the Fermi pockets, but the gap function is actually a 5x5 matrix. Its experimental implications 
O are also discussed. 
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\ The discovery of the high-Tc cuprate was epoch-making 
in that the material has opened a new avenue of the corre- 
lated electron systems and associated electron mechanism 
of superconductivity. However, we are still some way from 
a complete understanding, and a new class of materials has 
been desired, which would be a testbench for various ideas 
and concepts, both experimental and theoretical, in the 
electron correlation. So the discovery of superconductivity 
in the iron-based pnictide LaFeAsO doped with fluorine 
(LaFeAsOi_a;Fa^) discovered by Hosono's group[l] is par- 
ticularly welcome. Indeed, this is the the first non-copper 
family of compounds whose Tc exceeds 50 K in the subse- 
quent studies (Fig{T]). 

This immediately stimulates renewed interests in the 
electron mechanism of high Tc superconductivity. Espe- 
cially interesting, in our view, is that the way in which the 
correlation effects such as magnetism and superconductiv- 
ity appear in strongly correlated electron systems is very 
sensitive to the underlying band structure. We can envisage 
to exploit the "fermiology in correlated electron systems" 
to explore the possibility of manipulating superconductiv- 
ity by choosing the shape of the Fermi surface. [2j 

In this paper we first briefly look at unconventional su- 
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Fig. 1. (color online) Tc against calendar year for various classes of 
superconductors. 

perconductivity as a background. We then give a descrip- 
tion of the electronic structure and a possible mechanism 
for superconductivity for the iron-based compound. Theo- 
retically, the material turns out to have a curious multiband 
structure, with the Fermi surface comprising some pock- 
ets. In order to investigate the pairing mechanism, here we 
first construct an electronic model for the iron-based com- 



pound with the "downfolding" based on a first-principles 
electronic structure calculation. The minimal model turns 
out to contain all the five Fe d orbit als. We then apply the 
five-band random-phase approximation (RPA) to solve the 
Eliashberg equation. We conclude that a nesting between 
multiple Fermi surface (pockets) results in a development 
of peculiar spin fluctuation modes, which in turn realises an 
unconventional s (roughly extended s) pairing. [3j The re- 
sult is intriguing as one realisation of the idea of fermiology 
in general and the idea of the "disconnected Fermi surface" 
in particular [4J. Finally we shall summarise that, as op- 
posed the strongly correlated, half-filled, one-band cuprate 
as a doped Mott insulator, the iron compound is moder- 
ately correlated, dilute multi-band system. An outlook for 
the iron compound is also given. 

2. Unconventional superconductivity 

After the discovery of the cuprate, a number of other 
classes of superconductors have been discovered, among 
which are the ruthenate, the cobaltate, the Hf compound 
and the Ce compound. These have given a strong impe- 
tus towards the study of electron mechanism of supercon- 
ductivity. In the conventional mechanism, superconductiv- 
ity arises from the electron-electron attraction mediated 
by phonons, where Tc ~ 10 K is one order of magni- 
tude smaller than the phonon energy cod ^ 100 K. In 
the electron mechanism, superconductivity arises from the 
electron-electron repulsion, where spin and charge fluctua- 
tions mediate pairing with Tc ~ 100 K, two orders of mag- 
nitude smaller than the electronic energy Ep ^ 10000 K. 

The pairing from repulsive interactions itself is nothing 
strange, if we allow anisotropic pairing symmetries accom- 
panied by anisotropic gap functions, A{k). The point is, if 
we look at the BCS gap equation, 

AW = -i:n«.,<=')^ta„l.(|il). (1) 

k' 

where V{k^ k') is the pairing interaction, we can immedi- 
ately recognise that a repulsion acts as an attraction when 
A(/c) changes sign across the typical momentum transfer 
(usually dictated by the dominant spin fluctuation mode). 
This is what is happening in the cuprate with a d-wave 
pairing. 

The question, then, is "why is Tc ~ ~^^f in the elec- 
tron mechanism is so low?", as most transparently dis- 
played by Uemura's experimental plot for Tc against Tp [5j . 
Theoretically, there are good reasons why Tc is low: (a) 
the effective attraction mediated by the fluctuation is much 
smaller than the bare (repulsive) interaction, (b) Quasi- 
particles are short-lived due to large self-energy correc- 
tions from the electron correlation, and (c) the prerequisite 
anisotropy in pairing with the gap-function nodes, which 
usually intersect the Fermi surface, suppresses Tc- We can 
then think of how we can manipulate these factors. Sev- 
eral years ago Kuroki and Arita[4j have proposed that we 



can overcome the difficulty (c) by considering disconnected 
Fermi surfaces, on which we can pierce the nodes in between 
the Fermi pockets. Each pocket is then fully gapped, with 
opposite signs across the pockets. While the disconnected 
Fermi surfaces may first seem artificial, it has been recog- 
nised that they actually occur in real materials including 
TMTSF[6J and Co compounds [7J. 

As for the spatial dimensionality, Arita et al have 
shown that two-dimensional (layered) systems are more 
favourable than 3D ones for the spin-fluctuation mediated 
superconductivity for ordinary lattices [8j. This is because 
the fraction of the phase space over which the pairing in- 
teraction is appreciable is larger in 2D. The result, also 
obtained by Monthoux and Lonzarich[9j, agrees with the 
empirical fact that recently discovered superconductors 
(cuprates, Co compound, Hf compound, CeCoIns, etc) are 
mostly layer-structured. The iron-based compound has 
turned out to be no exception. 

3. LaFeAsO 

LaFeAsO has a tetragonal layered structure (Fig 12]), 
where each Fe layer comprises edge-shared FeAs4 tetrahe- 
dra with Fe atoms forming a square lattice. Since the bands 
around Ep have negligible LaO-layer components, we can 
concentrate on the Fe layer. The experimentally determined 
lattice constants are a — 4.03552Aand c = 8.7393A, with 
two internal coordinates z\,^ = 0.1415 and 2:as = 0.6512. 
In the experimentally obtained phase diagram in the orig- 
inal discovery [Ij and in an /iSR study [10jarXiv0806. 3533 
the (colinear; see below) SDW is taken over by the super- 
conducting phase as the F-doping level is increased. 




Fig. 2. (color online) Crystal structure of the LaFeAsO. The right 
panel picks up the FeAs layers. 

After the discovery of superconductivity, a surge of 
intensive studies on the material has ensued. Since it 
is impossible to survey them here, let us here men- 
tion only a few examples with an extremely partial list 
of references. In the iron-based family, various com- 
pounds have been shown to exhibit superconductivity, 
among which are NdFeAsOi_^F^[llJ. PrFeAsOi_^F^[12]. 



GdFeAsOi_^[13j, SmFeAsOi_^F^[14j, with Tc now ex- 
ceeding 56 K. Pressure effects have been studied for 
LaFeAsOi_^F^p^. NdFeAsOi_^[16j, etc. Superconduc- 
tivity has also been found in iron-based materials with 
different layered structures that include BaFe2As2[l7] and 
FeSeilSj. 

We can first itemise the experimental results that indi- 
cate unconventional SC: 

(i) NMR, where absence of coherence peak, 1/Ti ~ 

etc have been observed. [19j 

(ii) Electronic specific heat 7 ex etc has been 
observed. [20J 

(iii) Point-contact tunnelling, where zero-bias peak, in- 
dicative of a sign change across Fermi pockets, etc 
have been observed. [2 Ij 

(iv) /iSR is used to measure the superfluid density. [22j 
Especially interesting is the Uemura plot, in which 
the iron compound is shown to reside more or less on 
the Uemura line, so in this sense the iron compound 
is at least as good as the cuprate. 

(v) Raman spectra indicate a weak electron-phonon 
coupling. [23] 

(vi) i^c2, Hci indicate a two-gap behaviour. [24] 

(vii) Photoemission spectroscopy, [25j and more recently 
ARPES depict the Fermi surface, gap, etc[26j. 

Experimental results for magnetic properties, especially 
interesting in relation to the mechanism of superconduc- 
tivity, include the following: 

(i) Mossbauer studies indicate a magnetic transition 
temperature Tn ^ 140 K with a magnetic moment 
M = 0.25 - 0.35/iB (at T = 4 K) [27l28] . 

(ii) Neutron diffraction studies indicate a colinear 
(stripe) SDW (see Figl5|) in Fe moments with Tn ^ 
134 K, M - 0.35/iB (at T = 8 K)[29j. 

(iii) /iSR suggests a spontaneous muon spin precession 
below Tat = 138 K[30j. 

(iv) NMR shows Tn = 142 K[31j. 

(v) Transport measurements indicate an anomaly around 
T = 140- 150 K[32J. 

(vi) X-ray diffraction shows a structural phase transition 
(tetragonal orthorhombic) at 150 K[33j, which 
is close to, but somewhat deviates from, the magnetic 
transition temperature. 

Theoretically, if we first look at the periodic table, cop- 
per, in its Cu^+ state, has a 3dP electron configuration. 
This makes the cuprate one-band, nearly half-filled system, 
where we have a spin 1/2 at each site with strong anti- 
ferromagnetic fluctuations. By contrast, Fe, located in the 
middle of the table, has a 3d^ configuration, so multiple 
d bands tend to be involved. A point of interest is that 
we usually associate iron with ferromagnetism rather than 
with superconductivity [34], and how the multiband struc- 
ture for the iron compound gives rise to superconductivity. 

We then definitely require the electronic band structure 
to begin with. 

- Several groups have obtained the band structure [35'36'37] 

- Theoretical suggestiongs on magnetic properties include 



an LDA+DMFT study which indicates no local moment 
for the Hund's coupling Jh < 0.35 eV[38 . 

- An ab initio (cRPA) calculation estimats Jh — 0.2 
eV.[39j 

- An itinerant charachter of the compound suggested from 
several papers[41j. 

- Sensitivity of the electronic structure against the Fe-As 
distance has been suggested in a maximally localised 
Wannier study that incorporates As 4p orbitals[42j. 

4. Model construction 

First two papers that study the mechanism of supercon- 
ductivity are by Mazin et al[40j and by Kuroki et al[3]. The 
former starts from the Fermi surface and construct a an ef- 
fective RPA treatment for the spin susceptibility xo (q) / [1 ~ 
J{q)Xo{q)] to predict the colinear SDW and an extended 
s-wave pairing. The latter starts with a model building via 
downfolding to actually solve Eliashberg's equation for the 
full five-band model. Here we briefly describe this approach. 
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Fig. 3. (color online) Left panel: The dispersion of the three bands 
that intersect Ep in LaFeAsOi-xF^ with 10 % doping (x = 0.1). 
Right panel: Fermi surface (with the inter-layer hopping ignored). 

The approach starts from a band structure (FigslSl HI) 
with the density- functional approximation with plane- wave 
basis [43 for constructing the model. For that, we have 
employed a standard procedure with maximally localized 
Wannier functions (MLWFs)[44j. These MLWFs, centred 
at the two Fe sites in the unit cell, have five orbital sym- 
metries (orbital 1: d^z'^-R'^^ 2: dxz^ 3: dyz^ 4: dx'^-v'^^ 5: 
dxy)- Dne to the tetrahedral coordination of As, there are 
two Fe atoms per unit cell, and X, F, Z refer to those in the 
original unit cell. The two Wannier orbitals in each unit 
cell are equivalent in that each Fe atom has the same local 
arrangement of other atoms. We can thus take a unit cell 




Fig. 4. (color online) The band structure of the five-band model in 
the unfolded BZ, where the inter-layer hoppings are included. Orbital 
character is indicated for the bands and for the Fermi surface (right 
panel, with the original (dashed lines) and the unfolded (solid) BZ 
indicated). Bottom panels depict (two equivalent realisations of) the 
colinear SDW. 

that contains only one orbital per symmetry by unfolding 
the Brillouin zone (BZ), and we end up with an effective 
five-band model on a square lattice, where x and y axes are 
rotated by 45 degrees from X-Y ^ to which we refer for all 
the wave vectors hereafter. We define the band filling n as 
the number of electrons/number of sites (e.g., n = 10 for 
full filling). The doping level x in LaFeAsOi-a^F^^ is related 
to the band filling as n = 6 + x. 

The five bands are heavily entangled as shown in FigHJ 
which refiects hybridisation of the five orbitals due to 
the tetrahedral coordination of As Ap orbitals around Fe. 
Hence we conclude that the minimal electronic model re- 
quires all the five bands. [35] In FigHJ the Fermi surface for 
n = 6.1 (corresponding to x = 0.1) obtained by ignoring 
the inter-layer hoppings is shown in the two-dimensional 
unfolded BZ. The shape of the Fermi surface is now being 
experimentally detected with ARPES.[26j 

The Fermi surface consists of four pieces (pockets in 
2D when the inter-layer hoppings are neglected): two con- 
centric hole pockets (denoted as ai, 0^2) centred around 
{kxi ky) = (0, 0), two electron pockets around (tt, 0) {/3i) or 
(0,7r) (/32), respectively, ai {/3i) corresponds to the Fermi 
surface around the FZ (MA) line (in the original BZ) in 
Singh's band calculation. [36j In addition, we notice that a 
portion of the band near (tt, tt) is fiat and touches the Ef 
at n = 6.1, so that the portion acts as a "quasi Fermi sur- 
face (7)" around (7r,7r). The orbital character is again en- 
tangled: a and portions of [3 near BZ edge has mainly dxz 
and dYz character, while the portions of /3 away from the 
BZ edge and 7 have mainly dx'^-Y'^ orbital character. 

5. Five-band RPA results for magnetism and 
pairing 

Having constructed the model, we move on to the RPA 
calculation in the 2D model. For the many-body part of 
the Hamiltonian, we consider the standard interaction 




Fig. 5. (color online) 5-band RPA result for the spin susceptibility 
for 10% doped LaFeAsOi-a^Fo; for U = 1.2, U' = 0.9, J = J' = 0.15, 
and T = 0.02 (in eV). Bottom panels depict the colinear SDW for 
k = (7r,0), (0,7r). 

terms that comprise the intra-orbital Coulomb [/, the 
inter-orbital Coulomb U' ^ the Hund's coupling J, and 
the pair- hopping J^ The modification of the band struc- 
ture due to the self-energy correction is not taken into 
account. For the five-band model. Green's function is a 
5x5 matrix in the orbital representation, while each of 
the spin and orbital susceptibilities Xii,i2,i2,,iA{h = 1, •••,5) 
is a 25 X 25 matrix. The Green's function and the effective 
pairing interactions, obtained from the susceptibilities, are 
plugged into the linearised Eliashberg equation, and the 
gap function as a 5 x 5 matrix is obtained along with the 
eigenvalue A. 32 x 32 /c-point meshes and 1024 Matsubara 
frequencies are taken. We find that the spin fiuctuations 
dominate over orbital fiuctuations as far as /7 > U' ^ so we 
can characterise the system with the spin susceptibility. 
We denote the largest eigenvalue of the spin susceptibility 
matrix for iujn = as Xs{k). The gap function matrix at 
the lowest Matsubara frequency is transformed into the 
band representation by a unitary transformation, and its 
diagonal element for band i is denoted as (j)i{k). 

Let us look at the result for the spin susceptibility Xs for 
U = 1.2, U' = 0.9, J = J' = 0.15, andT = 0.02 (ah in units 
of eV) in Fig|5l The spin susceptibility has peaks around 
{kx^ky) = (tt, 0), (0, tt) (i.e., a colinear SDW), which refiects 
the Fermi surface nesting around ~ (7r,0), (0, tt) across a 
and P and also P and 7 in FigHl The spin structure is con- 
sistent with the neutron scattering experiments in which 
the colinear structure is observed at low temperatures for 
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Fig. 6. (color online) Bottom panels: Diagonal elements of the 
gap function matrix in the orbital representation for the XZ^ Y Z 
(left) and - y2 (right) in the 5-band RPA for 10% doped 
LaFeAsOi-xFcc for the same parameter set as in the previous figure. 
Nesting vectors are indicated in the top panels. 

the undoped case.[29j 

For superconductivity, we show in Fig l6] the diagonal el- 
ements of the gap function matrix for orbit als (iyz, dzx^ 
and (ix2-y2. The gap is a curious extended 5-wave,[40 3^ 
where the gap changes sign across the a — or — 7 nest- 
ing vector ~ (tt, 0), (0, tt) at which the spin fluctuations 
develop. The sign change is analogous to those in models 
studied by Bulut et a/.,[46j and Kuroki and Arita[4j. To be 
more precise for the multiband system, the magnitude of 
gap of the (ix2-y2 orbital turns out to be large compared 
to other orbitals, which indicates that the dx'^-Y'^ orbitals 
play the main role in the superconductivity. The extended 
s (~ coskx + co^ky in k space) should look in real space as 
depicted in FiglTl 

However, we have to be careful in analysing the gap, 
since we have a multiband system at hand, where the band 
character changes along the Fermi surface (FigH]). When 
we concentrate on the main component {dx'^-Y'^ orbitals) 
of the gap function, the magnitude of the gap varies along 
the (3 Fermi surface, which is because the character of the 
Fermi surface changes, as we approach the BZ boundary, to 
dvz^ dzx for which the gap is smaller. More importantly, 
we should always be aware that the gap function is a 5 x 5 
matrix, where off-diagonal elements are rather significant 
due to the heavy entanglement of the bands. This should 
have experimental implications for various quantities. An 
importance of matrix character of the gap is illustrated in 
FiglHlthat displays the quantity {(j)(j)^)4^4, where (j) is the gap 
matrix and 44 denotes the diagonal element of band 4. As 
shown in the figure, this quantity is fully gapped over the 
entire BZ. A remnant of the nodal lines of the diagonal el- 
ement appears as a dip that intersects the P Fermi surface, 
which may look like a node e.g. at not low enough temper- 
atures. The degree of the variation of the gap may be de- 
termined experimentally from e.g. tunnelling spectroscopy 



Cuprate Present study 




Fig. 7. (color online) Gap function symmetry is schematically shown 
for the cuprate (left panels) and the present result for the iron 
compound (right), in the /c-space (top panels) and in the real space 
(bottom). 
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Fig. 8. (color online) 5-band RPA result for the gap function, 
t)44 for the same parameter set as in the previous figures. 

or ARPES. Also interesting is how the penetration depth, 
A ~ ns/m* with Us being the superfluid density in the one- 
band system, should be given in the multiband system. 

6. Conclusion 

To summarise, we have very briefly reviewed the iron 
compound, and then described a theory in which a flve- 
band electronic model is constructed as a minimum micro- 
scopic model. Applying a flve-band RPA to this model, we 
have found that spin fluctuation modes around (tt, 0), (0, tt) 
develop due to the nesting between disconnected Fermi sur- 
faces. Based on the linearised Eliashberg equation, we have 
concluded that spin fluctuation modes realise an uncon- 
ventional, extended s-wave pairing, where the diagonal el- 
ements of the gap matrix change sign across the nesting 
vector. The d^s-ys orbital is the main component in the 
superconductivity gap function, although the gap function 
as a 5 X 5 matrix may be important. 

Overall, while the cuprate is strongly correlated {U > 
W with W the band width), one-band and nearly half- 



filled system, the iron compound is moderately correlated 
{U ~ iy)[39j, multi-band and dilutely filled system. This 
poses a challenging problem of whether the iron compound 
can exceed the cuprate. Relevance to experimental results, 
especially /iSR including the Uemura plot, is also a pressing 
future problem. 

The work described here is a collaboration with Kazuhiko 
Kuroki, Seiichiro Onari, Ryotaro Arita, Hidetomo Usui, 
Yukio Tanaka and Hiroshi Kontani. Numerical calculations 
were performed at the facilities of the Information Technol- 
ogy Center, University of Tokyo, and also at the Supercom- 
puter Center, ISSP, University of Tokyo. This study has 
been supported by Grants-in-Aid for Scientific Research 
from MEXT of Japan and from the Japan Society for the 
Promotion of Science. 
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